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Abstract 

A comprehensive parameterization of the colour singlet heavy quark free energy above T c is 
given, using the lattice data in quenched (Nf = 0) and unquenched (Nf = 2 and Nf = 3) QCD. 
The corresponding (temperature dependent) potentials thus obtained are then inserted into the 
Schrodinger equation for the charmonium and the bottomonium in the deconfined phase of QCD. 
The solution of the equation provides an estimate of the melting temperature and of the radii for 
the different cc and bb bound states. 
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I. INTRODUCTION 



The anomalous suppression of the J/if> production in heavy ion collisions, which has 
been experimentally observed in the depletion of the dilepton multiplicity in the 

region of invariant mass corresponding to the J /if} meson, was proposed long time ago as 
a possibly unambiguous signal of the onset of deconfinement Indeed in Ref. jjj it is 
argued that charmonium states can only be produced in the first instants after the nucleus- 
nucleus collision, before the formation of a thermalized QGP. Then, in their path through 
the deconfined medium, the original cc bound states tend to melt, since the binding (colour) 
Coulomb potential is screened by the large number of colour charges. This, in turn, produces 
an anomalous (with respect to normal nuclear absorption) drop in the J /if} yields. 

In this picture it is implicitly assumed that, once the charmonium dissociates, the heavy 
quarks hadronize by combining with light quarks only (recombination leading to a secondary 
J /i/j production is neglected). This assumption is certainly justified at the SPS conditions, 
due to the very small number of cc pairs produced per collision (N c5 ~ 0.2 in a central 
collision), but at RHIC (iV c5 ~ 10) and LHC (N c5 ~ 200) energies it is no longer warranted 

I 

Moreover in a hadronic collisions only about 60% of the observed J/if)'s are directly 
produced, the remaining stemming from the decays of excited charmonium states (notably 
the Xc and the if)'). Since each cc bound state dissociates at a different temperature, a model 
of sequential suppression was developed, with the aim of reproducing the J/ if) suppression 
pattern as a function of the energy density reached in the heavy ion collision (the highest 
temperatures and energy densities being reached in the most central collisions) p, (], B, S, E| . 
SPS experimental data for Pb-Pb collisions at different centralities seem indeed to support 
the dissociation pattern predicted by this model l[Q|. 

Alternative mechanisms for the J /if) production, like the Statistical Coalescence Model 

. In the SCM one assumes again that all the cc pairs 
are produced in hard processes at the initial stage of the collision. Any heavy quark bound 
state, if present, is assumed to melt in the QGP phase and the number of cc pairs in the 
fireball is considered fixed. At the chemical freeze-out open and hidden charmed hadrons 
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are then produced with multiplicity ratios 1 fixed by their masses, according to the laws of 
statistical mechanics. Hence, in such a scheme, the measured J ftp multiplicity is not related 
to the presence of charmonium bound states in the plasma phase, but to the statistical 
hadronization of the initially produced c(c) ( ant i-) quarks. This may lead, at LHC energies, 
to a completely different picture characterized by an enhanced charmonium production even 
if all the cc bound states dissociate during the plasma phase. 

In any case the hypothesis that all the primary produced J/ip's melt during the QGP 
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lifetime is hardly realized at SPS conditions. Hence models have been developed 
attempting to account both for the initial state production, eventually subject to in-medium 
dissociation, and for the thermal production at the hadronization. 

Concerning the heavy q nar k onia in the QGP phase, recent lattice data Q U Q in 
quenched approximation (hence neglecting effects arising from virtual processes involving 
dynamical fermions), which display narrow peaks for the charmonium spectral functions in 
the pseudoscalar and vector channels (even up to T ~ 2T C ), seem to point to the existence of 
heavy quark bound states up to temperatures above T c 2 . Clearly these results, if confirmed, 
would entail striking experimental consequences. 

Actually, the meson spectral functions cannot be measured directly on the lattice. From 
the numerical simulations one gets the current-current correlation function along the (imag- 
inary) temporal direction on a finite number of points. Such a correlator corresponds to 



the convolution of the meson spectral function with a thermal kerne 



The spectral function 
a procedure called 



can then be obtained only indirectly. With this aim in Refs. Q, Q, Q 
Maximum Entropy Method (MEM) has been adopted. Clearly, an independent check of 
the results obtained with MEM appears desirable. This indeed is our scope in the present 
paper. 

For this purpose, we first extract from lattice data a heavy quark potential accounting 



for thermal effects and then we solve t 



re Schrodinger equation for the charmonium (and 



bottomonium). As shown in Refs. 



y 0, y y 



l20l |. from the Polyakov loop correlation 



function it is possible to extract the free energy (in the different color channels) of a heavy 
quark- ant iquark pair placed at a distance r in a thermal bath of gluons and light dynamical 



1 Due to the low rate of inelastic reactions full chemical equilibrium cannot be reached by charmed hadrons: 
their total multiplicity measured at SPS stays well above the thermal value. 

2 Actually, in Ref. |b| the analysis has been done for ss mesons. 
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fermions. Once a good parameterization of the color singlet free energy is obtained, the 
entropy and internal energy contributions can be disentangled. Since the quarks acting as 
static sources of the color field are considered infinitely heavy, the internal energy coincides 
with the potential. The latter is then inserted into the Schrodinger equation, from which 
the binding energy of the different stable states — if there are any — and their evolution 
with the temperature are obtained. 

Indeed a clear distinction between the QQ free and potential energies is necessary in order 
to get a reliable estimate of the quarkonium dissociation temperature T d in the different spin- 
parity channels. 



In Refs. U m XM 



221 ]. where the colour singlet free energy was directly inserted into the 



Schrodinger equation, the dissociation temperatures T d = 1.10T C ?1 and T d = 0.99T C 
were found for the J/ip, all the other charmonium states melting well below T c . 
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On the other hand, in Ref. |23(, where a parameterization of the lattice color singlet 
potential (in the quenched approximation) was used, the temperature T d for the spontaneous 
dissociation of the J/ip was estimated to occur at about 2T r; a value even larger than the one 
obtained from the spectral analysis performed in Refs. ^,14, 3]. Also results for different 
charmonium and bottomonium states have been reported in Ref. j^]. The case Nf = 2 was 



addressed in Ref. 



241 ] . where the J/ip meson was found to be bound till T ~ 2.7T C . 



Actually, even if the potential supports the existence of bound states, other physical 
processes may lead to the dissociation of the quarkonium. First, if the QQ binding energy is 
lower than the temperature — and assuming that the quarkonia have reached the thermal 
equilibrium with the plasma — a certain fraction of their total number will be thermally 
excited to resonant states according to a Bose-Einstein distribution: such a process is referred 
to as thermal dissociation j^lL l^. Furthermore, the collisions with the gluons and the light 
quarks of the plasma may lead to the collisional dissociation of the quarkonium. In this 
connection, the reaction g + J/tp — > c + c was studied in detail in Ref. j^J. Hence, in 
spite of the presence of a bound state solution of the Schrodinger equation till T ~ 2T C , the 
J/ip turned out to be really stable only up to temperatures lower than 2t c 2^. Clearly the 
two processes above mentioned are not encoded in the Polyakov loop correlation function, 
where the heavy quarks act as unthermalized static sources of the color field, and have to 
be accounted for a posteriori. 

Of course, if the process g + J/ip — > c + c can lead to the dissociation of the charmonium, 
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the same reaction can also occur in the opposite direction. Hence a consistent calculation of 
J/ip multiplicity implies the solution of a kinetic rate equation integrated over the lifetime 
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of the QGP phase in which both processes (dissociation and recombination) enter 
To carry out this detailed balance calculation the knowledge of the J / ip binding energy and 
wave function in the thermal bath turns out to be an important input. This is of relevance 
because, as mentioned above, the usual assumption in considering the J/tj) suppression as 
a signature of deconfinement is that its production can occur only in the very initial stage 
of the collision. Really, if at SPS the role played by recombination is numerically negligible, 
this is no longer true at RHIC as pointed out in Ref. 

In any case in this paper we limit ourselves to check the existence of bound state solutions 
of the Schrodinger equation. A quantitative study of the dissociation and recombination 
processes is left for future work. 

Here we take advantage of all the available lattice data, obtained not only in quenched 
QCD (Nf = 0), but also including two and, more recently, three light flavors. We are then 
in a position to study also the flavor dependence of the dissociation process, a perspective 
not yet achieved by the parallel studies of the spectral functions, which are, as already 
mentioned, only available in quenched QCD. 

The present paper is organized as follows. In Sec. |H] we present a parameterization of 
ihe color singlet free energy lattice data for the cases Nf = jl7|, Nf = 2 [2(| and Nf = 3 
2fi| , from which the heavy quark potential is obtained. In Sec. IIHI we solve numerically 
the associated Schrodinger equation at different temperatures for the charmonium and bot- 
tomonium states, thus determining their dissociation temperature. Finally, in Sec. IIVI we 
present our conclusions. 



II. PARAMETERIZATION OF THE LATTICE DATA 



In this section we provide a unified parameterization of the lattice data for the color 



singlet QQ free energy iq in the case of quenched 



17], 



2-flavor 20] and 3-flavor |2J QCD. 



The lattice findings are shown (in dimensionless units) in Fig. ^ 

For what concerns the critical temperature we assume the values T r = 270 MeV (Nf = 0), 
T c = 202 MeV (N f = 2) and T c = 193 MeV (N f = 3) given in Ref. I20J. 

The free energy on the lattice is defined up to an additive normalization constant, which 
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FIG. 1: A sample of the available data for F\(r, T)/\fa at different values of T/T c as a function of 

n 

the separation r\fa of the QQ sources. Data are taken from Ref. [17j (Nf = 0), Ref. 20j (Nf = 2) 
and Ref. |26J (Nf = 3). Dimensionless units are employed, with y/a = 420 MeV (Nf = and 
N f = 2) and y/a = 460 MeV (JV/ = 3). 



has to be fixed using some physical constraint. In Refs. it has been normalized 

to match, at the shortest available distance for each temperature, the T = heavy quark 
potential. Such a normalization amounts to make the (reasonable) assumption that thermal 
effects become negligible at very short distances. 
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In Refs. [iij, [20j the zero temperature heavy quark potential has been determined through 
a best fit procedure of the available T = lattice data with a Cornell-like parameterization 3 : 

V(r) 4 a r - 

a representing the string tension. The values a = 0.195(1) for the case Nf = and a = 
0.212(3) for the case Nf = 2 are given in Ref. [Q, where the value \/a = 420 MeV is 
employed to translate the lattice results into physical units. 

In Ref. a similar parameterization (Cornell potential plus a 1/r 2 term to mimic the 
effects of asymptotic freedom at the shortest distances reachable on the lattice) is employed 
for the case Nf = 3. Note that in that work the free energy is provided directly in physical 
units: however, for the sake of comparison, we show it in dimensionless units, using for the 
string tension the value \fa = 460 MeV extracted from the parameterization of the T = 
potential given in Ref. |26|. 

In past calculations [6. l22j ] the free energy has been often identified with the 

heavy quark potential and inserted directly into the Schrodinger equation. However, a 
better candidate for a more appropriate finite temperature potential is given by the internal 
energy of the QQ system, defined by the well known relation 

F = U-TS, (2) 

where 

V = (3) 

is the internal energy and 

OF 

is the entropy. One can see from the data in Fig. ^ that the role played by the entropy is 
more relevant at large distances. 

Getting the internal energy from the free energy involves a derivative of the latter with 
respect to the temperature: it is thus clear that one needs an accurate parameterization of 
the temperature dependence of F. 



3 Note that the 1/r term accounts for two different physical processes: the perturbative one-gluon-exchange 
at short distances and the transverse string fluctuations at large distances. 
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In order to establish a suitable form for this parameterization of the lattice data, we first 
consider the two limits in which the underlying physics is supposed to be known. 

At very short distances (r <C ^-/T) thermal effects are negligible and the colour singlet 
free energy is dominated by the perturbative one-gluon exchange with the typical behaviour: 

(5) 

the coupling a depending only upon the QQ separation. 

On the other hand, for T 3> T c , the large distances (r ^> 1/T) behaviour of the free energy 
is expected to be dominated by the exchange of a resummed electrostatic gluon leading to 
the expression 27]: 

F 1 {r,T) ~ -l^le- mD ^ + F 1 (r = oo,r). (6) 

In this limit, the coupling a is a function of the temperature and mr,{T) is the Debye 
screening mass arising from the dressing of the electrostatic gluon. 

In the two above limits, the running of the coupling is determined by a Renormalization 
Group Equation (RGE), allowing to express a = g 2 /^ as a known (at least in the weak 
coupling regime) function of an energy scale //. In the short distance limit the relevant 
energy scale is given by the inverse of the distance (// ~ 1/r), while for large separations the 
major role in setting the scale is expected to be played by the temperature (/i ~ T). 

Actually, in order to solve the Schrodinger equation for heavy quarkonia, one really needs 
a parameterization of the free energy covering the whole range of distances. For this purpose, 
on the basis of Eqs. (0) and © it appears convenient to cast the dependence of F\ on r and 
T into the following functional form 4 : 

3 r 

In order to recover from the above the short and large distance limits given by Eqs. (0) and 
(jBJ) one can, e.g., assume a to depend on the following combination of r and T: 

a(r,T) = ol{\l = c r /r + c t T), (8) 



4 Note that the mass M(T) appearing in this exponential will not necessarily coincide with the Debye 
screening mass of Eq. (0 , the latter being determined by fitting only the large distance data. 
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FIG. 2: The colour singlet free energy resulting from our fitting procedure compared to the lattice 
data at different temperatures above T c . 



where a(fi) is obtained by solving the RGE, while c r and q are numerical coefficients to be 
fixed through a best fit of the data. 

If supported by the data, Eq. (|SJ) would allow to interpolate between the short distance 
regime (r <C 1/T), where fi ~ 1/r, and the long distance one (r ^> 1/T), where on the 
contrary /i ~ T. 

We now describe in detail the fitting procedure we have employed. The data for the 
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colour singlet free energy, taken from Ref. jl7| (Nf = 0), Ref. Q (Nf = 2) and Ref. 2(| 
(Nf = 3), are displayed in Fig. [T]in dimensionless units, namely y = F\j ^fa and x = r^/a, 
for different values of the temperature both below and above T c . For each value of T > T c 
the data have been parameterized with the following (four parameter) fitting function: 

4^S_ e -a 3 x + ao with ~ <^l +a2 (9) 

3 x x 

where for we use the RGE result obtained with the two-loop QCD beta-function quoted 
in Appendix El the dimensionless variable jl being identified with the ratio /x/Aqcd- 

The fitting procedure yields a very mild dependence on T for the parameters a\ and 
Actually, on the basis of the above discussion, one would have expected the coefficient a 2 
to scale linearly with the temperature, but our finding might just signal that the range of 
temperatures spanned here is still not in the asymptotic regime T 3> T c . On the other hand, 
a constant value for a x is in agreement with the ansatz of Eq. (jSJ). This has suggested the 
following procedure: a weighted average of the values obtained for a\ has been performed, 
yielding a x = 0.2719(2) for N f = 0, a x = 0.2687(7) for N f = 2 and a x = 0.2354(17) for 
Nf = 3. We have then fitted again the data keeping ai fixed and using only a , a 2 and (Z3 
as free parameters. This parameterization works well — yielding values of y 2 per degree of 
freedom of the order of 1 at all the temperatures — and it is compared to the data in Fig. |21 

The finite temperature lattice data are limited to distances r > 0.1 -v- 0.2 fm. Since 
the data examined in this work have been normalized assuming that at short distances 
thermal effects are negligible, we should check that our parameterization does not introduce 
any sizable (and spurious) temperature dependence for small values of r, remaining in this 
region close to the T = perturbative potential (we remind that at T = free energy 
and internal energy coincide). In this respect our choice of refitting the data keeping the 
coefficient ai fixed fulfills this requirement. 

The behaviour of our parameterization of the free energy for r < 0.1 fm is displayed, 
at three different temperatures, in Fig. El where we indeed see that spurious short distance 
thermal effects appears negligible. It is gratifying (and somewhat surprising) that our curves 
seems to interpolate smoothly between the T = Cornell potential and the short distance 
perturbative potential. The Cornell curve reported in the figure is the one employed in 
Re, B to Ex the notation of toe free ene rgy at the «e re at te-np— A S 
already discussed at the beginning of this section it was obtained by fitting zero temperature 
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FIG. 3: The short distance behaviour of the colour singlet free energy resulting from our fits for 
different values of T/T c . The dependence on the temperature turns out to be very small at such 
short distances. We also display for comparison the T = potentials given by different schemes: the 
potentials employed to normalize the data at short distance 20, 2f|, the potential obtained in the 
Nambu string model and the one and two-loop short distance perturbative potentials 



28, 



29 



3- 



lattice data which cover distances r > 0.05 fm. Of course for shorter distances the Cornell 
parameterization cannot account for running coupling effects (asymptotic freedom) and the 



perturbative calculation 



28. 2i 



should provide more reliable results. More details on 



the one- and two-loop perturbative potential used in Fig. |3] are given in Appendix |XJ 
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The curves V^ tr i ng reported in Fig. prefer to the QQ potential 

V(r) = -~+ar (10) 

obtained in the Nambu string model, the term 1/r arising from the quantum fluctuations of 
the flux tube in the transverse directions. 

As a next step we have provided a parameterization of the T-dependence of do, a>2 and 
03, which allows one to perform analytically the derivative of F\ with respect to the tem- 
perature (actually, for the calculation of the QQ binding energies the knowledge of ao is 
unnecessary: we show also this contribution to the potential for completeness). Lacking 
compelling physical hints to the form of the T-dependence, we have sought for the sim- 
plest and smoothest expressions yielding a satisfactory interpolation of the values calculated 
from the fit to the lattice data (see Appendix El for details). The parameters and their 
interpolations are displayed in Fig. HJ 

The value of Oo(T) at each temperature is essentially fixed by the normalization of the 
data, which has been determined using similar procedures in all the lattice calculations. 
Notably this parameter, once plotted as a function of T/T c , turns out to be fairly close for 
all values of Nf. On the other hand, the remaining two parameters display some quantitative 
flavor dependence, but within the same qualitative pattern. An exception to this behavior is 
represented by 02 (T) for Nf = 3, although one should note that the three-flavor case is the 
one where the parameterization © has the largest uncertainties. Moreover, the variation of 
ci2 (T) with T is actually magnified by the scale of the figure, being generally around 20% for 
the range of temperatures considered here. As we shall see the resulting effective potential 
is rather robust with respect to these variations. 

Note also that the strongest variation of the parameters with the temperature is confined 
around T very close to T c . Since extracting the internal energy from Fi involves a derivative 
with respect to T, this is the temperature domain where the sensitivity to the details of 
the parameterization might be high. For this reason we felt it safer to use the resulting 
potentials at temperatures larger than T c , say for T > 1.05T C , where they turn out to be 
stable with respect to changes in the parameterization. 

Finally, we recall that the quenched QCD data are actually available also at temper- 
atures much larger than those for the unquenched cases: here we limit our analysis up to 
about twice the critical temperature, since this is the range where lattice calculations with 
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FIG. 4: Temperature dependence of the parameters ao, a,i and 03 employed in fitting the color 
singlet free energy with Eq. <P|). The solid lines are interpolating functions. 



two and three flavors are available and, moreover, since this is the range of interest for the 
problem of the J/ip dissociation. Here we just remark that an analysis of the Nf = data 
at larger temperatures (till T ~ 5T C ) with the parameterization of Eq. Q shows that the 
fitting parameters maintain, also at these high temperatures, the trend displayed in Fig. EJ 
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III. QQ BOUND STATES IN QUARK-GLUON PLASMA 



The QQ free energy obtained in lattice calculations has been used in the past as an 



input for the QQ potential energy in the non-relativistic Schrodinger equation 
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More recently |l7l . l23l . |2J] it has been recognized that the QQ potential energy can be more 



appropriately identified with the QQ internal energy (see Eq. (j2J)) 5 . 

Once the temperature dependence of the color-singlet free energy F\ has been param- 
eterized, the corresponding color-singlet internal energy U\ is easily obtained and one can 
define an effective potential 



V 1 (r,T) = U 1 (r,T)-U 1 (r^oo,T) 



(11) 



to be used in the Schrodinger equation 

V 2 



2/i 



+ Vx(r,T) 



1>(r,T) = e(T)i/>{r,T), 



(12) 



where \i is the reduced mass of the QQ system. 



5 Actually, the author of Ref. [2^] tries to disentangle, in the total internal energy Ui(r,T), the gluon and 
the QQ contributions. This leads to a lower value for the J/ip dissociation temperature with respect to 
the result ~ 2T C obtained with the full internal energy. 
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FIG. 6: Binding energy of charmonium states above the deconfinement temperature. The lines 
show the results for m c = 1.3 GeV; the grey areas around the lines display the variations for 
1.15 GeV< m c <1.35 GeV. 



In Fig. El we display, at two different temperatures, the effective QQ potentials that we 
have obtained by using the parameterizations discussed in the previous section. At high 
temperature the form of V\ is the one typical of a screened Coulomb potential, with the 
Nf = 3 potential providing a stronger attraction; at temperatures close to T c the shape of 
V\ appears somewhat distorted because of the strong temperature dependence of a 2 and a 3 
(see Fig. HJ), the Nf = 3 potential being still more attractive. Remarkably, the behaviour 
of the color singlet potential energies obtained with our procedure appears qualitatively in 
agreement with the one given in Refs. HQ, 

where the internal energy has been directly 
calculated on the lattice for Nf = and Nf = 2, respectively. 
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We now employ the effective potential previously derived for the study of the charmonium 
and bottomonium spectroscopy above the critical temperature. 

For what concerns the charmonium states we plot the values of the binding energies 
(Fig. EJ) and of the mean square radii (Fig. |7J) of the different bound state solutions of 
Eq. (|12j) as functions of T/T c . For the sake of comparison the results obtained for a different 
number of dynamical fermions (Nf = 0, 2, 3) are plotted in the same panel. 

In our analysis we have chosen m c = 1.3 GeV for the charm quark mass; however, we 
have also studied the sensitivity of our results to variations of m c in the range given in the 
Particle Data Group (PDG) listings j^, namely 1.15 GeV< m c < 1 .35 GeV. For reference, 
in Fig. [7| the arrows point to the T = mean square radii obtained in a potential model 
calculation j^j]. Note that potential models at T = typically follow a different philosophy: 
they fix the quark mass and the confining potential parameters to reproduce the mass of the 
lowest states, whereas here the potential is provided by lattice calculations and the quark 



mass is the running mass in the MS scheme. 

The IP and 2S states turn out to melt at temperatures 1.1T C < Td < 1.15T C . On the 
contrary the J /if) stays bound up to temperatures l.7T c < Td < 2.3T C , the precise limits 
depending upon Nf. The lower bound (T d ~ 1.7T C ) refers to the quenched case and appears 
in striking agreement with the limiting value obtained in Ref. ^| through the study of the 
J/ip and rj c spectral functions. 

Note that the free energy measured on the lattice provides a spin averaged result of 
the singlet and triplet channels: hence, in Figs. El and U\ the J / ip and t] c mesons appear 
degenerate. This degeneracy is expected to be removed by a short range spin-spin force, 



whose effect, at T = 0, is often treated perturbatively assuming a contact interaction 



3: 



8vra s (/i) 

H ss = — <T • (TqS[r). (13) 

9 m q mq 

Again, the short range nature of this force makes plausible the assumption that it is not 
affected by thermal effects and that it can be employed also at finite temperatures. Then, 
for the J/ij) and r\ c energy shifts one gets: 

A£, w = ^^%(0, T)| 2 (14a) 
9 m r m r . 



AE Vc = -^^%(0,T)| 2 , (14b) 
9 m c m e 

where for a s (/i) we have used the two- loop expression of Eq. (|A5j) . evaluated at [A — m c . As 
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FIG. 7: Mean square radii of the charmonium states above the deconfinement temperature. The 
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one can see, the temperature dependance of this term stems entirely from the value of the 
cc wave function at the origin. 

We display in Fig. |S] the binding energies of the J/ if) and r] c mesons above the decon- 
fmement temperature. The spin-spin contribution gives rise to a mild reduction of the 
dissociation temperature, without altering the qualitative features of the results. For in- 
stance, for Nf = 2 at T = 1.05T C one gets a J/ip — r] c splitting of 145 MeV (the experimental 
value at T = is about 117 MeV), whereas at the J/ip melting temperature, T = 1.8T C , 
one gets a splitting of 13 MeV (although, strictly speaking, in this case the perturbative 
expansion is no longer justified). 

In Figs. El and El we plot the binding energies and the mean square radii of the bottomo- 
nium bound states above T c . We have chosen m;, = 4.3 GeV for the bottom quark mass, but 
again we have studied the sensitivity to variations of rrib in the range of the PDG listings, 
4.1 GeV< rrih<4A GeV. The arrows in Fig.^lpoint to the results of T = potential model 
calculations |33j]. Of course, for heavier quark masses more states than in the charmonium 
case survive above the deconfmement temperature. 

Both for the charmonium and for the bottomonium states the binding energy gets smaller 
and the mean square radius gets larger as the temperature grows. This should reflect into 
a huge increase of the elastic cross section a(QQ — > QQ) at low energy. In fact, for a 
small relative momentum k of the two particles, when in the partial wave expansion the 
contribution of the s-wave is dominant, the presence of a / = state near zero binding 
energy leads to a dramatic increase of the cross section as k — > 0, according to the formula: 

47T 

ai =° k z fTW (15) 

where [i is reduced mass of the system and e the energy of the state near zero binding, 
no matter whether it is positive (virtual level) or negative (bound state). If something 
analogous happened also for the heavy-light states (for which lattice data are not available 
yet) this would clearly accelerate the thermalization of the heavy quarks. 

From Figs. EJ El an d El the effective potential obtained with Nf = 3 appears more attrac- 
tive. Indeed the critical temperature T c in the three cases Nf = 0, 2, 3 assumes different 
values, hence a naive comparison might not be so meaningful. Furthermore the finite tem- 
perature Nf = 3 lattice data available so far are affected by larger errors, so that our 
parameterization presents larger uncertainties with respect to the other two cases. 
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FIG. 9: Binding energy of the bottomonium states above the deconfinement temperature. The 
lines show the results for m& = 4.3 GeV; the grey areas around the lines display the variations for 
4.1 GeV< m b <4.4 GeV. 
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IV. CONCLUSIONS 



The picture of the deconfined phase of QCD for values of the temperature slightly ex- 
ceeding T c — that is the ones accessible in the heavy-ion experiments presently performed 
at RHIC — has substantially evolved in the last few years. 

For large values of the temperature (T> 3T C ) a description of the QGP in terms of a 
gas of weakly interacting quasiparticles |35( appears reliable (even in a regime where the 
coupling g is not so weak) and supported by the lattice data for the QGP thermodynamics 

On the other hand, for temperatures up to T ~ 2T C (namely the ones currently accessible 
in the experiments), the matter resulting from the heavy-ion collisions is nowadays often 



described 



24j as a strongly interacting QGP (sQGP). In particular, a striking feature of the 



QGP matter obtained at RHIC is its hydrodynamical behaviour |37|, |38( (characterized by 
a very low viscosity), which manifests itself in particular in the elliptic flow observed in 
non-central collisions: the plasma obtained at RHIC seems to behave as a nearly ideal fluid 
whose expansion is driven by pressure gradients. 

In order to explain the very small mean free path of the plasma particles required by the 
hydrodynamical scenario, a picture of the matter obtained at RHIC in terms of a system of 
lundreds of loosely bound states of quasiparticles (qq, qg, gg...) has been recently proposed 



24j . In such a framework one has to resort to some assumptions. Gluons and light quarks 



are treated as quasiparticles endowed with quite heavy thermal masses obtained from lattice 
calculations and the potential felt by them in the different color channels is got from the 
lattice QQ free energy under the hypothesis of Casimir scaling. 

The presence of such a pattern of bound states in the range of temperatures T c < T < 
2T C has been recently questioned on the basis of the analysis of the correlation between 
baryon number and strangeness j^j]. The evaluation of the correlation coefficient Cbs — 
—3(BS)/(S 2 ) should allow to discriminate between a scenario in which the relevant degrees 
of freedom in the QGP are weakly interacting quark and gluon quasiparticles or loosely 
bound states of the latters. The lattice data on the off-diagonal quark-number susceptibilities 
available so far seem to favour the first hypothesis. 

Here we followed a different approach, starting from the case of a heavy QQ pair placed in 
a thermalized QGP and extracting their interaction from the available lattice calculations. 
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Indeed, we have exploited the lattice data for the heavy quark free energies to get information 
on the existence of cc and bb bound states above the deconfmement transition. We have 
examined the cases Nf = [l^l, Nf = 2 (^(J and also Nf = 3 26 1, where lattice data 
are getting available. For the color singlet free energy we have adopted a parameterization 
which accounts for the effects of asymptotic freedom at short distances and displays an 
exponential screening at large distances. From the free energy we have then extracted the 
heavy quark potential, to be inserted into the Schrodinger equation. The latter has been 
solved numerically for the two interesting cases of charmonium and bottomonium. 

For what concerns the charmonium, we have found a dissociation temperature Td ~ 2T C 
for the IS states (r] c and J/ij)). Also the excited states IP and 2S appear to melt, but at 
temperatures slightly exceeding T c . 

On the other hand, the bottomonium spectrum displays a much larger number of bound 
states above T c . In particular its ground state turns out to remain bound in the whole range 
of temperatures covered by our parameterization. By extrapolating the latter at larger T 
we get a dissociation temperature Td ~ 4 6T C , depending upon the number of dynamical 
fermions. At fixed T/T c , both for the charmonium and for the bottomonium, the system 
turns out to be more bound when the number of light flavors is increased. 

Remarkably, in the range of temperatures of experimental interest covered in this paper 
a number of loosely bound QQ states exists. Since the existence of states near zero binding 
energy entails a huge increase of the elastic cross sections at low relative momenta, this can 
help the approach to thermal equilibrium also of the heavy quarks. 

An analogous situation occurs in the cross-over from the BCS theory of superconductivity 
to the Bose-Einstein condensation, much explored lately in ensembles of alkaline fermionic 
atoms confined in a magnetic trap. Indeed, by smoothly changing the external magnetic 
field B, one induces a smooth change of the energy of the two interacting atoms: if this 
energy is c 



occurs [41, 



ose to the one of a Feshbach resonance, a huge increase in the cross-section 
42, 13, 44}]. The striking analogy between B and T as control parameters is 



self-imposing. 

In connection with the previous discussion on heavy quark thermalization, it would be 
highly desirable to have lattice data also for the Qq states, which are not yet available. 

A phenomenological model has been recently proposed assuming the existence of res- 
onant (not bound) D- and B- meson states above T c (4^. In this scheme the transverse 
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momentum distributions of the charmed quarks (anti-quarks) would approach their ther- 
mal equilibrium value much faster, due to isotropic resonant scattering on light anti-quarks 
(quarks). Indeed, recent PHENIX results for azimuthally averaged transverse momentum 
spectra of single electrons arising, in Au-Au collisions, from the decay of D- and B-mesons 
seem to be compatible with a thermalization scenario 4|| . 

Furthermore, preliminary results from PHENIX and STAR on the elliptic flow of the 
above electrons tend to support a picture in which, due to strong rescattering, charmed 
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50j . Hence, 



quarks reach thermal equilibrium and follow the flow of the fireball 
an extension of the present approach to the Qq (qQ) states might shed light on the mechanism 
of thermalization of the heavy quarks. 

Coming back to the present results, they offer the relevant possibility of evaluating the 
charmonium (bottomonium) multiplicity produced in heavy ion collisions at the experimen- 
tal conditions of SPS, RHIC and LHC 

Indeed, a reliable estimate of the dissociation temperature of the different quarkonia is 
essential in order to predict how many of them, after being produced in the hard initial 
processes, survive in the QGP phase thus contributing to the final measured yields. The 
latter of course also contain, e. g., the charmonia that might be thermally produced during 
the hadronization process by cc recombination. 

In particular, the knowledge of the in-medium quarkonium wave function and binding 
energy would allow a hopefully reliable estimate of the gluon-dissociation cross section of 
the J /if) (and hence, via a detailed balance analysis, also of the cross section for the inverse 
process). Thus a kinetic rate equation accounting for both the dissociation and recombina- 
tion processes can be tackled: its solution, depending on how many cc pairs are produced 
in the initial state, should provide the number of J/ip present at the end of the QGP phase. 
These issues will be addressed in future work. 
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APPENDIX A: RUNNING COUPLING AND PERTURB ATIVE POTENTIALS 
AT SHORT DISTANCES 



In this ap pen dix we collect the formulas resulting from the perturbative calculations of 



Refs. 



28, 



291 . |30j of the heavy quark potential at short distances, where asymptotic freedom 
guarantees that the perturbative approach is justified. 
The static QCD potential turns out to be given by: 
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In the above r' = rexp(7#) (je is the Euler-Mascheroni constant), C A = 3, C F = 4/3, 
T F — 1/2 and (3 = Cr(3) (the Riemann function of argument 3); moreover, ot^ip) is the 



QCD running coupling in the MS renormalization scheme coming from the solution of the 
RGE: 



da s (fi 2 



n=0 
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n+1 



(A4) 



a s (/j, 2 ) dlog/j, 2 

The one-loop perturbative potential is by obtained keeping only the first two terms of 
Eq. (|A1|) and employing for cxj^g(fx) the result 
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which arises from the solution of Eq. (|A4J) with the two-loop beta function (i.e. considering 
only the first two terms of the series). Indeed the coefficients (3q and (3\ do not depend on 
the renormalization scheme and are given by: 
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(A6) 
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On the other hand, for the two-loop perturbative potential one has to keep all the terms up 
to order displayed in Eq. (jAl|) . The evaluation of the running coupling should include 
also the three-loop coefficient $2 of the beta-function, which is no more renormalization 
scheme independent. In evaluating the two-loop curve in Fig. |3]we have used Eq. (13) of 
Ref. 0, 

where the MS scheme is employed. 
Clearly, Eq. ()A1|) still depends on the parameters [i and Aqcd- The scale /i ~ r _1 is quite 
arbitrary: we chose /i = [r exp^^)] -1 as usually done in the literature. For what concerns 
Aqcd we employed the value suggested by a recent lattice collaboration jsi| Aq CD = 261 
MeV for all the three cases (N f = 0, 2, 3). 



APPENDIX B: T-DEPENDENCE OF THE FIT PARAMETERS 

In this appendix we show the functions employed in fitting the temperature dependence 
of the parameters ao(T), a2(T) and 03 (T) entering into the functional form we have adopted 
for the QQ free energy (see Sect. ITTj) . 

As already mentioned in the text, since we have no phenomeno logical or theoretical hints 
to the functional form of the T-dependence of these parameters, we have tried to use the 
simplest expressions yielding a smooth fit and a "reasonable" \ 2 ( sa Y) °f the order of a few 
units). 

For ao(T) and CI2 {T) we have used the following forms: 

4 0) ^ 0) exp[-4 0) x] 

a ° = — — ' (B1) 

a - A ( 2 ) 4 2) +4 2) * 2 + * 4 (B2) 
A\ I} + Afx + x 3 

where x = T/T c and the Af^s and A\ 2 ^s are fit parameters. 

Only in the case of a 3 (T) we have used an expression inspired by perturbative QCD, 
namely the one for the Debye mass (although, rigorously this parameter does not represent 
the Debye mass): 

( n f \ 1/2 / (3) 4 3) 4 3) 4 3) 4 3) \ 

a3= { 1 + -f) X [ A » 1 + ^ + ^ + ir + ^J^ioo P (*), (B3) 

where g^-ioop/^ 71 = a Jis ^ s the 2-loop coupling constant obtained by replacing /x/Aqcd with 
4.8826a; ^ (27tT c /Aq C d)x, having used T c = 202 MeV and A QCD = 261 MeV. In the case 
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Nf = 3, because of the larger errors affecting the parameterization of the free energy, we 

(3) 

made the fit of a% stiffer by setting A\ = 0. 
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